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1.  INTRODUCTION 


Coherent  laser  detection  and  ranging  (LADAR)  systems  are  currently  being 
developed  for  the  acquisition  and  tracking  of  extended  targets.  Although 
heterodyne  detection  of  the  optical  signal  generally  provides  significant 
signal-to-noise  ratio  enhancement  over  incoherent  detection  for  Infrared 
wavelengths,  the  average  signal  (and  resultant  slgnal-to-nolse  ratio)  may  be 
significantly  degraded  by  phase-front  distortions  resulting  from  target 
speckle  and  atmospheric  turbulence  effects,  unless  proper  system-design 
precautions  are  observed.  Similarly,  monopulse  tracking  performance  of  LADAR 
systems  can  be  degraded  unless  phase-distortion  effects  are  understood  and 
compensated.  The  objective  of  this  study  was  to  develop  analyses  and  computer 
codes  that  permit  reliable  predictions  of  phase-distortion  effects  on 
heterodyne  LADAR  performance  for  realistic  system  parameters,  target 
configurations,  and  engagement  scenarios. 

This  final  report  summarizes  the  results  obtained  from  15  July  1980 
through  15  July  1983  concerning  the  efficiency  of  the  heterodyne-detection 
process  resulting  solely  from  extra-instrumental  phase  effects  (i.e., 
unavoidable  instrument  inefficiencies  are  neglected).  An  extended  period  of 
performance  was  provided  to  permit  the  comparison  of  analytical  predictions, 
developed  under  this  contract,  with  experimental  results  that  were  to  be 
provided  under  other  programs.  The  experimental  results  developed  during  this 
time  period  were  of  limited  value,  so  the  emphasis  of  this  study  was  re¬ 
directed  toward  (a)  refining  the  analyses  previously  developed  (MDC  Final 
Report  Q0714,  Laser  Radar  Analyses),  (b)  developing  and  documenting  a  general 
heterodyne-efficiency  prediction  code,  and  (c)  providing  illustrative 
calculations  that  can  be  referenced  for  future  experimental  and  system  design 
activities.  This  report  summarizes  the  results  of  these  efforts. 


‘1.  SUMMARY  OF  RESULTS 


Significant  results  of  the  contract  study  for  the  period  IS  July  1080 
through  15  July  1983  are  summarized  below.  Detailed  analytical  results  are 
provided  in  Section  5,  and  Illustrative  numerical  results  are  given  In 
Section  6. 

•  A  general,  flexible,  computer  code  was  developed  to  compute  the 
laser-radar  heterodyne-efficiency  degradation  that  results  from 
phase-front  distortions  introduced  by  both  atmospheric  turbulence 
and  target-dependent  speckle  effects.  The  code  permits  speckle 
calculations  of  either  a  Schell-model  target  or  a  general  target 
having  an  arbitrary  projected  cross-section  distribution. 
Computations  can  be  performed  assuming  either  a  circular-aperture  or 
a  Cassegrain  receiver. 

•  A  comparison  of  computed  heterodyne-efficiency  factors  with 
experimental  data  obtained  at  the  General  Electric  Co.  (G.E.) 
Electronics  Laboratory  indicated  general  agreement  in  trends  and 
reasonable  agreement  in  magnitude  for  a  diffuse,  scattering  target 
but  poor  agreement  In  magnitude  for  a  specular  target.  Uncertainty 
In  experimental  conditions  render  this  comparison  Inconclusive . 

•  Good  correspc  ence  of  the.  predictions  of  the  computer  code  was 
demonstrated  with  published  results  of  a  more  restrictive  study  that 
analyzed  only  turbulence  effects  coupled  with  restrictive  receiver 
parameters . 

•  Illustrative  results  were  computed  for  the  speckle  degradation  in 
heterodyne  efficiency  that  would  be  expected  for  both  a  simulated 
tactical  aircraft  and  a  Schell-model  target,  as  a  function  of  a 
range  and  beam  radius,  to  emulate  the  detection  efficiency  that 
could  be  expected  in  an  operational  environment. 

•  A  test  panel  configuration  that  could  be  employed  for  future 
experimental  tests  was  developed,  and  efficiency  computations  were 
performed  on  the  test  array  that  demonstrate  the  feasibility  of 
experimental  tests  of  predictions  concerning  coherence  effects  and  a 
phase-quadrature  component. 
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3.  CONCLUSIONS 


Salient  conclusions  resulting  from  this  study  during  the  reporting  period 
are  summarized  below. 

•  The  heterodyne-efficiency  computer  code  that  was  developed  is 
sufficient  for  systems  analyses  and  planned  experimental  tests 
concerning  the  effects  of  atmospheric  turbulence  and  target  speckle 
on  the  heterodyne  efficiency  of  laser  radars. 

•  The  high  confidence  in  the  validity  of  the  model  used  for  heterodyne 
efficiency,  and  the  resultant  computer  code,  indicates  that  the  lack 
of  agreement  with  G.E,  experimental  data  probably  r« suited  from 
inadequately  characterized  experimental  effe;;s  (such  as  beam 
jitter)  and  inadequate  data  acquisition  and  processing  techniques. 

•  The  rapid  decrease  in  the  heterodyne  efficiency  with  decreasing 
range  in  the  near-field  will  probably  dictate  either  the  use  of 
focusing  optics  (in  either  transmitter  or  receiver)  or  corrective 
data  processing  to  maintain  a  slgnal-to-noise  ratio  of  a  laser  radar 
system  equivalent  to  that  at  intermediate  ranges. 

•  Predicted  negative  values  of  the  laser-radar  heterodyne  efficiency 
correspond  to  a  detected  signal  that  is  out  of  phase  with  the 
intermediate  (heterodyne)  frequency. 

•  The  actual  signal-to-noise  ratio  of  a  heterodyne  laser-radar  system 
can  be  increased  by  decreasing  the  transmitted  beam  divergence,  but 
there  exists  a  maximum  benefit  that  is  dependent  on  the  particular 
transceiver  configuration. 
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4.  RECOMMENDATIONS 


•  Laser  radar  experiments  should  be  conducted  using  the  test  panel 
configuration,  suggested  In  Section  6.4,  to  verify  the  effects  of  speckle 
on  heterodyne  efficiency  and  the  correspondence  of  predicted  negative 
efficiencies  with  a  phase-quadrature  component. 

•  The  analyses  developed  In  this  Investigation  should  be  extended  to  permit 
the  study  of  speckle  effects  on  monopulse  tracking. 

•  The  results  of  this  study,  which  yield  the  effective  carrier-to-nolse 
ratio,  should  be  Integrated  Into  a  more  general  analysis  of  the  complete 
slgnal-to-noise  ratio,  that  determines  the  detection  probability, 
required  for  acquisition  studies. 
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5.  ANALYSES 


To  provide  a  convenient  reference  for  the  equations  employed  in  the 
generalized  laser-radar  efficiency  code  (section  7),  relevant  formulae  derived 
In  the  previous  investigation  (Reference  1)  are  summarized  below  In  Section 
5.1.  A  corrected  form  of  the  limits  of  Integration  that  arise  In  the  Incohe¬ 
rent  form  of  the  heterodyne-efficiency  Integral  Is  provided  In  this  summary, 
as  well  as  some  comments  concerning  the  Interpretation  of  the  computed 
values.  Section  5.2  provides  a  brief  derivation  of  the  modification  of  the 
heterodyne-efficiency  formulae  that  result  when  the  heterodyne  efficiency  of  a 
receiver  having  Cassegrain  optics  Is  desired. 

5.1  Heterodyne  Efficiency  Equetlons 

Beginning  with  the  equations  for  the  ensemble-averaged  photodetector 
current  correlation  function  of  a  heterodyne  LADAR  system  developed  in 
Reference  2,  the  previous  Investigation  showed  that  the  maximum  power  slgnal- 
t.o-nolse  ratio  heterodyne  efficiency  factor  Is  given  by 

r («')  »*„(!> 

V  jj  ■ 

‘Sjp 

where  y  Is  the  spatial  degree-of-coherence  function  at  the  receiving  aperture 
plane,  u^b  is  the  normalized  spatial  variation  of  the  back-propagated  local- 
oscillator  wave,  Ar  is  the  receiving  aperture  area,  and  A^p  is  the  back- 
propagated  aperture  area  (i.e.,  the  minimum  areal  substance  of  either  the 
receiving  aperture  or  the  back-propagated  local-oscillator  wave,  whichever  is 
less) . 

To  permit  specific  calculations  of  LADAR  efficiency  factors,  the 
generalized  transceiver  system  depicted  in  Figure  1  is  considered,  having  the 
following  parameters: 
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Figure  1.  Transceiver  parameters  and  coordinates. 


A  -  local-oscillator  aperture  radius, 

P  »  entrance  pupil  radius, 

f^  -  focal  length  of  local  oscillator  lens/mirror  (lens  3),  and 
fr  «  focal  length  of  receiving  lens/mirror  (lens  2). 

These  parameters  were  used  in  Reference  1  to  derive  an  excellent  approximation 
for  the  back-propagated  local-oscillator  field  using  the  definitions 
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The  derived  expression  for  the  back.'  propagated  local-oscillator  field 
facilitates  the  evaluation  of  the  heterodyne  detection  efficiency, 

Equation  (1),  subject  to  various  assumptions  concerning  the  mutual  intensity 
functional  form  of  the  received  optical  signal.  If  the  incident  field 
correlations  can  be  approximated  by  a  Schell-model  function  of  the  form 


2  -|<a|!|2-Kx*jr|2) 

U  e  e 


(5) 


where 


(6) 


rc  is  the  coherence  radius,  L/^  is  the  amplitude  radius,  and  is  the  average 
phase-front  radius  of  curvature,  it  is  straightforward  to  show  that 
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Assuming  an  on-axis  Schell-model  target,  the  incident  field  correlation 
parameters  appearing  in  Equations  (5—8)  can  be  calculated  from  the  following 
expressions  (Reference  3) 
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where  a',  appearing  in  Equations  (10-12),  is  the  Equation  (6)  parameter  that 
results  from  the  target  amplitude  and  phase-front  radii  (L^  and  ),  R  is 
the  target  range,  and  pQ  and  are  the  atmospheric  and  target  coherence 
lengths,  respectively. 

A  more  general  expression  for  the  heterodyne  detection  efficiency  is 
obtained  by  relaxing  the  constraint  that  the  target  radiance  is  Gaussian 
(yielding  the  Equation  (5)  functional  form  for  the  mutual  intensity).  If  the 
target  has  a  single  glint  point  and  an  arbitrary  diffuse  reflectivity 
distribution,  the  degree-of-coherence  function  will  be  given  by* 
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where 
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range  to  target, 

transverse  oordinate  to  line-of-sight  at  the  target  plane, 


d2o(p) 


differential  projected  scattering  cross-section  per  unit  area, 


aj  and  &2  are  the  Gaussian  radii  of  curvature  at  the  glint  point,  A  is  a 
factor  that  relates  the  relative  strength  of  the  coherent  (glint)  and 
incoherent  components,  and  the  denominator  of  Equation  (13)  results  from  the 
normalization  requirement  of  the  degree-of-coherence  function. 

The  parameter  p0  in  Equation  (13)  is  the  turbulent-atmosphere  spherical- 
wave  coherence  length  £pQ  =»  [3/8  x  1.455  k^RC^ ]“3/5 >  where  is  the 
refractive  index  structure  constant^  describing  the  effect  of  the  assumed 
homogenecus  atmospheric  turbulence  on  the  propagation  or  partially  coherent 
light.4 

Recognizing  that  the  generalized  degree-of-coherence  function  of  Equation 
(13)  consists  of  the  sum  of  a  coherent  and  incoherent  component,  the 
calculated  heterodyne-efficiency  factor,  given  by  Equation  (1),  likewise  can 
be  split  into  the  sum  of  two  components.  The  maximum  heterodyne  efficiency 
can  be  expressed  as 
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where  otar  is  the  total  target  incoherent  LRCS,  i.e.. 
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and  the  coherent  and  incoherent  efficiencies  are  determined  from  the  relation 
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the  integral  factors  QCoh/incoh  are  defined  by  the  integral  numerator  of 
Equation  (1).  Because  the  coherent  component  of  the  degree-of-coherence 
function  is  in  the  form  of  a  Schell-model  MIF,  the  coherent  component  of  the 
total  efficiency  can  be  expressed  in  the  form  of  Equation  (7)  with  the 
variable  aQ  replaced  by  the  variable  <*1  defined  as 
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The  incoherent  integral  factor  can  be  written  in  the  form 
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and  -  pr  +  2p~2, 

by  transforming  the  integration  limits  of  Equation  (14)  to  the  sum  and 
difference  coordinates  and  utilizing  the  convolution  Integral  expression  for 
the  Fourier  transform  of  the  product  of  two  functions.  The  notations 
pr  and  p^  employed  In  the  above  expressions  denote  the  real  and  imaginary 
parts  of  p,  respectively.  The  function  F( )  defined  by  Equation  (20) 
results  from  the  Fourier  transform  of  the  turbulent-atmosphere  spherical-wave 
MIF  truncated  by  the  finite  receiver  aperture  limits  and  Gaussian  local- 
oscillator  field  distribution. ^  Physically,  F  represents  the  spatial 
filtering  of  the  target's  spatial  power  spectrum  by  the  turbulent  atmosphere 
and  finite  limits  of  the  receiver.  The  Equation  (20)  expression  for  the 
filtering  function,  F,  differs  from  that  reported  in  Reference  1  (Equation 
68).  The  corrected  expression,  Equation  (20),  reflects  the  fact  that  the  sum 
and  difference  coordinates 

t  -  i  d  +  V) 

and 

t  -  t  -  V  (22) 

can  ba  expressed  In  terms  of  the  aperture  coordinates  %  and  %'  to  yield  a 
restriction  on  their  combined  limits,  viz., 

2f*2  +  I  f 2  -  t2  +  t'2  <2  Rffl2.  (23) 

Thus,  the  difference  coordinate,  that  is  the  integration  variable  in  the  inner 
integral  of  Equation  (18),  is  restricted  to  the  range 

t2  <  4  (R2  -  t2),  04) 

which  is  reflected  in  the  parameters  appearing  in  the  filter  function  in 
Equation  (20). 
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To  numerically  evaluate  Equation  (18),  it  is  desirable  to  employ  an 
approximate  analytical  form  for  the  filtering  function  F.  It  is  demonstrated 
in  Reference  l  that  an  excellent  approximation  to  Equation  (20)  results  from 
the  function 


F(S') 


a'  > 

U'  < 


(25) 


(26) 


(27) 


The  parameter  P  appearing  in  Equations  (25)-(27)  is  the  integration  limit  of 
the  difference  coordinate,  i.e., 


(28) 


which  replaces  the  back-propagated  aperture  radius,  in  the  corresponding 
expression  (Equations  (79)— (72))  of  Reference  1.  It  should  also  be  noted  that 
Equation  (18)  differs  from  Equation  (67)  of  Reference  1  by  a  factor  of  1/2. 
This  factor  results  from  the  observation  that  because 


(29) 


when  is  transformed  to  -  £~,  the  apercure  coordinates  of  Equation  (29)  are 
interchanged.  Thus,  a  circularly  symmetric  integration  over  the  variable 
in  Equation  (18)  yields  a  double  mapping  over  the  coordinates  %  and  . 


It  Is  easy  to  demonstrate,  using  the  approximate  forms  for  the  filter 
function,  that  the  incoherent  detection  efficiency  approaches  unity  as  R^j  goes 
to  zero  and  approaches  zero  in  a  limiting  fashion  as  goes  to  infinity,  as 
required  from  the  physics  of  the  problem.  Computed  results  presented  in 
Section  6.2  demonstrate  the  correspondence  of  the  general  expression, 

Equation  (18),  with  the  more  restrictive  formula,  Equation  (7),  for  the 
special  case  of  a  Schell-model  source.  The  resultant  confidence  in  the 
approximations  employed  for  the  spatial  filtering  function,  F,  coupled  vith 
the  comparisons  of  the  preceding  analyses'*  with  independent  investigations  of 
a  more  limited  scope®,  support  the  validity  of  Equations  (18),  (19),  and 
Equations  (25)-(27).  These  formulae  provide  a  prescription  for  calculating 
the  heterodyne  efficiency  of  detecting  a  target  of  rather  arbitrary  laser- 
reflectivity  properties  in  the  presence  of  atmospheric  turbulence.  Although 
other  investigators  have  performed  extensive  analyses  of  the  signal-to-noise 
degradation  of  laser-radar  imagery  resulting  from  speckle  effects  and 
atmospheric  turbulence^,  a  similar  investigation  of  heterodyne  efficiency  has 
not  been  heretofore  reported. 

The  user  of  the  above  formulae  should  note  several  points  concerning 
their  derivation  and  interpretation.  First,  the  derivation  assumes  that  the 
mutual  intensity  function  (which  Is  a  statistical  measure  derived  from  an 
ensemble  average)  adequately  describes  the  time-averaged  phase-front 
correlations  that  give  rise  to  the  time-averaged  laser-radar  signal.  This 
assumption  is  predicated  on  a  form  of  the  ergodic  hypothesis®  which  holds  when 
the  phase  deviations  produced  by  the  scattering  are  large  and  the  scattered 
field  is  quasihomogeneous  so  that  a  spatial  average  samples  the  stochastic 
variations  of  the  system.  In  this  event,  the  spatial  variations  of  the 
scattering  geometry  that  change  with  time,  encountered  for  most  laser-radar 
systems,  yield  time-averaged  signals  that  should  be  in  good  accord  with  the 
prediction  of  an  ensemble-averaged  statistic.  Second,  the  derivation  is 
predicated  on  the  frequency  integrated  signal  power  which  is  given  by  the 
current  variance.*-  Thus,  the  derived  relations  describe  the  heterodyne 
efficiency  for  the  signal  power  included  in  all  spectral  components  of  the 
laser-radar  pulse-form  when  the  pulse  is  "on".  Finally,  it  should  be  noted 
that  the  numerator  of  the  derived  expression  for  the  heterodyne  signal-power 
efficiency,  Equation  (1),  is  such  that  negative  values  are  possible.  The  fact 
that  negative  values  of  the  mutual  intensity  function  can  dominate  the  inte- 
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grated  contribution  is  evident  from  the  Equation  (13)  form  of  the  degree-of- 
coherence  function,  which  shows  that  the  incoherent  contribution  can  be 
presented  in  the  form  of  a  Fourier  transform  of  the  differential  projected 
scattering  cross-section.  Clearly,  negative  values  of  the  Fourier  transform 
are  permissable  and,  in  fact,  this  eventuality  is  demonstrated  in  a  numerical 
example  provided  in  Section  6.  The  calculated  negative  values,  however,  pose 
a  problem  concerning  their  physical  interpretation  because  the  calculated 
quantity  is  a  power  ratio  which  should  always  be  positive.  Because  the 
derivation^  assumes  that  the  heterodyne  current  is  at  the  heterodyne 
Intermediate  frequency  (i.e.,  in  phase  with  the  difference  frequency  between 
the  laser  and  local  oscillator),  it  seems  reasonable  to  assume  that  negative 
values  correspond  to  the  out-of-phase  (quadrature)  component.  This  interpre¬ 
tation  is  supported  by  the  reasonable  physical  predictions  obtained  assuming  a 
physically  meaningful  power  value  for  negative  values  presented  in  Figure  14 
in  Section  6.  The  verification  of  this  physical  interpretation  of  the 
negative  efficiencies  as  a  quadrature  component  must,  however,  await  further 
experiments . 

5.2  Cassegrain  Receiver  Geometry 

The  equations  summarised  in  the  preceding  section  were  developed  for  the 
laser-radar  system  depicted  in  Figure  1,  where  an  unobstructed,  circularly 
symmetric,  receiver  collects  the  laser  light  scattered  from  the  target  through 
the  turbulent  atmosphere.  Because  the  phase-front  distortion  effects  on  the 
heterodyne  efficiency  are  calculated  from  the  phase  interference  with  the 
back-propagated  local-oscillator  wave  at  the  receiver  aperture,  the  foregoing 
equations  must  be  modified  to  treat  the  case  of  a  Cassegrain  receiver  where  a 
central  area  of  the  receiving  aperture  is  obscured.  The  following  development 
treats  the  necessary  modifications  for  the  Cassegrain  receiver  that  are 
implemented  in  the  general  heterodyne-efficiency  code  listed  in  Section  7. 

Because  all  of  the  heterodyne-efficiency  equations  are  of  integral  form, 
since  they  are  derived  from  the  basic  integral  equation  given  in  Equation  (1), 
it  is  sufficient  to  consider  a  general  integrand  of  two  spatial  variables, 
f(r,  r'),  integrated  between  limits  and  Rq  representing  the  inner 
obscuration  radius  and  the  outer  aperture  radius,  respectively.  It  Is 
straightforward  to  show  that 
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Thus,  the  integral  expression  for  the  Schell-modal  target  heterodyne 
efficiency,  Equation  (7),  and  the  coherent  component  heterodyne  efficiency 
(Equation  (17)  with  Equation  (7))  can  be  calculated  for  a  Cassegrain  receiver 
geometry  using  the  Equation  (30)  prescription  for  the  limits  of  integration 
appropriate  for  the  receiver.  It  is  also  straightforward  to  show  that  the 
efficiency  integral  is  normalized  by  a  factor  given  by 


ecoh/incoh 


(31) 


instead  of  Equation  (16)  in  the  case  of  Cassegrain  receiver.  The  notation 
employed  in  Equation  (31)  is  a  subscript  form  of  the  parameters  defined 
earlier,  viz.,  PQ  and  are  the  outer  aperture  radiuB  and  inner  aperture 
radius,  respectively,  and  Rffll  and  R^  are  the  respective  inner  and  outer  back- 
propagated  local-oscillator  beam  radii  (i.e.,  the  extension  of  Equation  (9)  to 
the  Cassegrain  geometry). 

Although  Equations  (30)  and  (31)  suffice  to  adapt  the  previous  formulae 
for  Schell-model  targets  and  the  coherent  component  of  scattering,  the  limits 
of  integration  of  the  incoherent  component,  Equation  (18),  must  be  further 
modified  because  this  equation  is  given  in  terms  of  the  sum  and  difference 
coordinates.  It  is  clear  from  the  above  development  that  the  limits  of 
Integration  on  three  Integrals  must  be  determined  in  terms  of  the  sum  and 
difference  coordinates  for  the  incoherent  scattering  contribution  for 
Cassegrain  geometries,  however  the  development  provided  below  demonstrates 
that  the  consideration  of  the  mixed  limits  integral  (the  last  integral  in 
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Equation  (30))  suffices  to  provide  the  limits  on  the  previous  two  integrals. 
The  mixed  limits  integral  of  two  spatial  variables  can  be  written 


EV 

J  J  d2U2V  f  {L  t') 


00 

ff  d2U2t  f(C,  CO  ©  (r J  -  C2)o  (rJ  -  C'2) 

—CD 


(32) 


in  terms  of  the  aperture  coordinates  and  the  Heaviside  functions,  G.  Because 
of  the  relations  (Equation  (22)),  between  the  aperture  coordinates  and  the  sum 
and  difference  coordinates,  it  follows  that 

if2  +  |  r2  -  t2  +  t'2  <  R2  +  R2  (33) 

for  the  mixed  limits  of  integration  case,  so  that 

t2  <  4(R2  -  r2),  (34) 


where 

R2  -  }  (R2  +  R2)  •  (35) 

In  a  similar  manner,  the  defining  equations  may  be  used  to  show  that 

t2  -  y  (I2  +  t'2  +  2%‘t'),  (36) 


t2  <  j  (R2  +  P2  +  2Rq  Rj  -  R+2,  (37) 

where 

R+  “  I  (Ro  +  Ri)*  (38) 

Therefore,  the  mixed  integration  limits  integral  may  be  written 
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x  e|A(R2  -  r2)  -  t2]  (39) 

in  terms  of  the  sum  and  difference  coordinates.  Equations  (35),  (38),  and 
(39)  show  that  the  incoherent  efficiency  scattering  integral  reduces  to  the 
form  previously  derived  for  an  unobscured  aperture  when  the  outer  and  inner 
integration  limits  are  equal.  Thus,  the  first  two  Integrations  reduce  to  the 
case  previously  treated,  Equations  (18)  and  (20),  with  replaced  by  and 
*ao  respectively.  For  the  case  of  mixed  integration  limits,  Equations  (18) 
and  (20)  become 


ff  d2r  d2r  i(f\  n  ®  (*+2  -  t2) 


Woh  ■  T  -  r2)  ff  d2?'  -  o] 
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(40) 


and 


u{r2-  f2  +'2/2p  }  r  - =2l  *  „  -"2/28 

F(^)  .  - ~ - e  11  J  du'M^"  2^RZ  -  r2Jlo(p>"/8)  e  * 


(41) 


respectively.  The  above  equations  provide  the  basis  for  the  modifications 
included  in  the  laser-radar  heterodyne-efficiency  code  of  Section  7  that 
permit  Cassegrain  receiver  efficiences  to  be  calculated  when  that  option  is 
selected. 
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6.  NUMERtCAL  RESULTS 


To  Illustrate  the  effects  of  atmospheric  turbulence  and  speckle  on 
heterodyne  detection-efficiency,  several  efficiency  functional  dependencies 
are  provided  below  that  were  numerically  computed  from  the  equations 
summarized  In  the  previous  section  as  implemented  in  the  generalized 
efficiency  code  presented  In  the  subsequent  section.  Specifically,  Section 
6.1  provides  a  comparison  of  the  calculated  results  with  the  only  applicable 
experimental  data  that  has  been  measured  to  date,  and  Sections  6.2  and  6.3 
provide  calculated  results  for  atmospheric  turbulence  degradations  of 
efficiency  for  a  clear  aperture  receiver  and  a  Cassegrain  receiver, 
respectively.  Section  6. A  provides  calculated  resells  of  the  speckle-related 
degradation  in  heterodyne  efficiency  for  a  conical  ’'arget,  a  representative 
tactical  aircraft  target,  and  a  test-panel  configuration  that  could  be 
employed  for  subsequent  tests. 

6.1  Comparison  of  Calculated  and  Experimental  Results 

Laser-radar  heterodyne-efficiency  factors  were  calculated  for  the  G.E. 
laser-radar  experimental  parameters  employed  on  13  ..  aly  1981  test.^  Assumed 
parameters  of  the  calculations  were  chosen  to  emulate  the  G.E.  experimental 
parameters  (i.e,,  local-oscillator  aperture  radlup  -  0.043  m,  Cassegrain 
obscuration  radius  »•  0,025  m,  Cassegrain  outer  radius  -  0.0381,  0.0509, 

0.0762,  0.0889,  and  0.C953  m,  receiver  focal  length  *  3.05  m,  local-oscillator 
focal  length  -  A. 06  m,  local-oscillator  beam  radlup  *  0,1  m,  target  range  =» 

609  ra,  and  transmit  beam  divergence  “0.6  rarad).  The  local-oscillator 
aperture  radius  was  chosen  to  emulate  the  finite  size  of  the  infrared  detector 
that  was  effectively  flooded  with  the  local-oscillator  beam  radiation.  A 
Schell-model  (Gaussian)  target  iEftjnsity  distribution  was  assumed  for  the 
diffuse  target  because  the  transmitted  Learn  was  prenumeably  a  TEMqq  mode  and 
the  target  was  flat  and  uniformly  diffuse.  The  calculated  results  are 
compared  with  the  experimental  data  in  Figure  2  as  a  function  of  the 
Cassegrain  aperture  diameter.  It  is  apparent  thac  the  calculated  results  for 
the  glint  target  exceed  measured  values,  while  the  diffuse  target  predictions 
are  somewhat  less  than  measured  results,  although  both  predictions  follow 
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Figure  2.  Calculated  and  measured  laser-radar  heterodyne  efficiency  factors  for  the  G.E.  laser-radar 
experimental  parameters  employed  on  the  13  July  1981  tests. 

observed  trends  (i.e»,  they  decrease  with  increasing  Cassegrain  aperture  size 
because  of  near-fleld  effects).  The  diffuse  target  results  would  be  closer  to 
observed  data  if  a  smaller  beam  divergence  were  assumed.  It  is  noteworthy 
that  Reference  9  emphasizes  that  the  beam  divergence  is  one  of  the  most  ill-' 
determined  parameters  of  the  experiment,  the  other  being  the  absolute 
heterodyne  efficiency.  Although  Figure  2  would  indicate  that  the  diffuse 
target  calculations  are  in  better  accord  with  the  experimental  data,  it  is 
apparent  that  systematic  errors  are  present  in  the  data  that  preclude  a  valid 
comparison.  These  errors  are  evident  by  the  fact  that  (1)  a  near-unity 
heterodyne  efficiency  was  not  measured  for  the  case  of  the  glint  target  with 
minimum  aperture  and  (2)  both  the  glint  and  diffuse  target  efficiencies  show 
an  anomalous  rise  for  the  5  in.  aperture  case.  Both  of  the  aforementioned 
experimental  anomalies  cannot  be  rationalized  with  the  known  physics  of  the 
detection  process.  It  is  possible  that  beam-jitter  problems  adversely 
affected  the  optical  alignment,  resulting  in  lower  measured  heterodyne 
efficiencies  for  the  glint  target.  The  diffuse  target  measurements  would  not 
be  as  subject  to  beam-jitter  problems  because  of  the  diffuse  target 
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scattering,  however,  the  speckle  effects  that  result  from  diffuse  scattering 
require  careful  data  recording  and  statistical  analyses  that  were  not 
adequately  performed  In  the  Reference  9  study.  Although  the  experimental  data 
generally  support  the  predicted  trends,  the  obvious  Introduction  of  artifacts 
renders  the  comparison  Inconclusive. 

6.2  Turbulence  Effects  on  a  Circular  Aperture  Heterodyne  Receiver 

Clifford  and  V.'andzurab  have  performed  calculations  of  the  effects  of  the 
turbulent  atmosphere  on  heterodyne  laser-radar  performance  assuming  a  focused 
beam  wave,  but  Ignoring  the  effects  of  a  finite  collection  aperture.  To 
demonstrate  the  correlation  of  the  predictions  of  this  analysis  with  the 
results  of  Clifford  and  Wandzura,  Figure  3  shows  computed  efficiencies  as  a 
function  of  integrated  turbulence,  assuming  the  unlikely  conditions  of 
Reference  6,  where  the  aperture  size  is  much  larger  than  the  Gaussian  beam 
radius.  The  integrated  turbulence  number  as  defined  by  Clifford  and  Wandzura 
is 
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Figure  3.  Calculated  heterodyne  efficiency  as  a  function  of  integrated  turbulence.  Calculations  assume 
A  =  P  =  5  m,  L0  =  0.3  m,  f?  =  fr  =  0.5  m,  and  X  =  10.6  pm. 
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where  the  factor  of  /I  results  from  the  different  definitions  of  the  Gaussian 
beam  radius  employed  in  their  investigation  and  this  one.  Figure  3  depicts 
the  functional  dependence  of  the  heterodyne  efficiency  on  the  integrated 
turbulence  for  three  different  range  values  corresponding  to  the  Fresnel 
numbers  (Q  ■  kL^/R)  of  approximately  0,  4,  and  16  used  in  the  Reference  6 
calculations.  The  Figure  3  results  were  computed  from  Equations  (7)  and  (17) 
and  agree  with  the  results  calculated  by  Clifford  and  Wandzura  using  the 
quadratic  atmospheric  structure-function  approximation^  assumed  in  this 
work.  However,  these  results  differ  slightly  from  results  computed  assuming 
the  more  exact  five-thirds  power-law  atmospheric  structure  function,  as 
pointed  out  in  Reference  6.  Larger  Fresnel  numbers  delay  the  decrease  of  the 
heterodyne  efficiency  with  increasing  turbulence,  although  the  initial 
efficiency  is  lower  because  of  near-field  effects  (i.e.,  poor  phase  matching 
on  the  photo-detector). 

Figure  4  shows  efficiency  calculations  similar  to  that  of  Figure  3; 
however,  the  receiver  aperture  was  assumed  to  be  the  same  size  as  the  local- 
oscillator  beam  radius  for  these  computations.  Figures  3a  and  3b  differ  only 
in  their  method  of  computation:  Figure  3a  was  computed  from  Equations  (7)  and 

(17) ,  whereas  Figure  3b  was  computed  from  the  more  general  Equations  (16), 

(18) ,  and  (20),  assuming  a  point-like  reflecting  object  having  a  radius  of 
only  1  mm.  The  discrepanci *s  between  Figures  4a  and  4b  result  solely  from  nu¬ 
merical  inaccuracies  in  computing  the  filter  function  of  Equation  (20). 
Numerical  computations  of  the  filter  function  become  increasingly  difficult 
for  larger  Fresnel  numbers  (Q),  and  the  approximations  of  Equations  (25)-(27) 
also  become  invalid.  The  Figure  4  results  are  particularly  interesting 
because  they  indicate  that  the  finite  receiving  aperture  has  the  effect  of 
further  delaying  the  decrease  of  the  efficiency  with  increasing  turbulence 
compared  with  the  beam-wave  local-oscillator  case  of  Figure  3,  and  a  small 
increase  of  heterodyne  efficiency  is  noted  for  certain  integrated  turbulence 
values  for  large  Fresnel  numbers.  Remarkably,  the  more-accurate  Figure  4a 
results  replicate  the  calculations  of  Clifford  and  Wandzura  for  the  beam-wave 
local-oscillator  configuration  using  the  more-accurate  assumption  of  a  five- 
thirds  power-law  atmospheric  structure  function,  although  a  quadratic 
structure  function  was  assumed  for  these  computations.  Clifford^  explains 
the  increase  in  heterodyne  efficiency  for  large  Fresnel  numbers  and  certain 
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Figure  4.  Calculated  heterodyne  efficiency  as  a  function  of  integrated  turbulence.  Calculations  assume 
A  =  P  =  Lq  =  0.3  m,  =  fr  =  0.5,  and  X  =  10.6  fim.  Figure  4a  is  computed  from  Eqs.  (7)  and 
(17)  and  Fig  4b  is  computed  irom  Eqs.  (18)  and  (25). 


turbulence  strengths  as  resulting  from  signal  increases  caused  by  the  focusing 
effect  of  large-scale  turbulent  eddies.  This  explanation,  however,  must  be 
re-examined  because  of  the  correspondence  of  Figure  4a  with  the  Reference  6 
results  and  the  differences  between  Figures  3  and  4. 

6.3  Turbulence  Effects  on  a  Cassegrain  Heterodyne  Receiver 

To  illustrate  the  effect  of  a  Cassegrain  receiver  geometry  (a  customary 
LADAR  design)  on  heterodyne  detection  efficiency,  Figure  5  shows  the  effici¬ 
ency  reduction  with  increasing  integrated  turbulence  for  the  same  conditions 
assumed  in  Figure  4  obscuration  having  a  radius  of  20  cm.  The  Figure  5 
results,  calculated  from  Equations  (7),  (17),  (30)  and  (31),  indicate  that  the 
central  obscuration  generally  has  the  effect  of  increasing  the  heterodyne 
efficiency  over  a  similarly  configured  unobscured  receiver,  although  anomalies 
are  apparent  for  the  high  Fresnel  number  calculation  (Q  *  16,  R  ■  1.1  km). 

This  conclusion  is  expected  because  of  the  larger  relative  size  of  the 
coherence  length  with  respect  to  the  unobscured  aperture  dimension.  If  the 
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Integrated  turbulence,  N0 

GPJI-llW-12 

Figure  5.  Calculated  heterodyne  efficiency  as  a  function  of  integrated  turbulence  for  a  Cassegrain 

receiver  geometry.  Calculations  assume  A  =  Lq  =  0.3  m,  fg  =  fr  =  0.5  m,  A  =  10.6  /xm,  and  an 
inner  (obscuration)  radius  ot  P,n  =  0.2  m  and  outer  (receiver  aperture)  radius 
of  pout  =  0-3  m• 

integrated  turbulence  number,  N0,  ic  redefined  to  more  accurately  reflect  the 
integration  area  of  the  Casaegt^  «.  receiver,  the  apparent  benefits  of  a 
Cassegrain  geometry  disappear. 


Speckle  Effects  on  Heterodyne  Efficiency 


Conical-Target  Scatter!™ 


Numerous  calculations  have  been  performed  using  the  analyses  summarised 
In  Section  5  to  assess  the  effects  of  phase— front  distortions  resulting  from 


target  speckle  on  the  heterodyne  efficiency.  Illustrative  of  the  types  of 
speckle  effects  that  are  encountered  are  the  cone  scattering  efficiency 
calculations  shown  in  Figure  6.  The  Figure  6  data  show  the  variation  in 


heterodyne  efficiency  as  a  function  of  the  aspect  angle  of  a  conical  target 
having  a  metallic  rough  surface  such  that  the  bidirectional  reflectivity  of 
the  surface  is  that  shown  in  Figure  7.  The  assumed  bidirectional  reflectivity 


function  is  typical  of  many  rough  surfaces  and  was  used  to  ca’culate  the 
projected  conical-target  cross  section  per  unit  area  (d  a/ dp  )  with  the  aid  of 
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Cone  backscatter  angle  (deg)  opJM199-13 

Figure  6.  Calculated  heterodyne-detection  efficiency  as  a  function  of  viewing  angle  for  a  conical  target 
having  a  base  radius  of  0.25  m,  a  height  of  1.0  m,  and  a  surface  finish  characterized  by  the 
bidirectional  reflectivity  distribution  depicted  in  Fig.  7.  The  fixed  parameters  of  the  calculation 
are  A  =  P  =  2,5  cm,  fg  =  fr  =  25  cm,  Lq  =  5  cm,  C^  =  10~**  m~^,  and  range  (R)  values 
are  annotated.  11 


Figure  7.  Laser  backscatter  bidirectional  distribution  fuuciions  for  hypothetical  metallic  rough-surface 
materials  having  two  scales  of  roughness  characterized  by  the  ROSSCO  parameters;  ko  -  0.2, 
kH  =  5.0,  kg  =  3.0,  s  =  0.6,  m-1.0,  n  =  100.0,  and  cr  =  0. 
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the  ROSSCO  code  summarized  in  Reference  il.  The  heterodyne  efficiencies, 
plotted  in  Figure  6  for  various  ranges,  were  computed  from  Equations  (18), 

(25),  (26),  and  (27)  using  the  four-fold  numerical  integration  employed  in  the 
generalized  efficiency  code  of  Section  7  that  has  been  shown  to  be 
convergent.  The  Figure  6  data  differ  slightly  from  the  corresponding  results 
depicted  in  Figure  25  of  Reference  1  because  of  the  modified  integration 
limits  employed  in  the  more  recent  calculations,  although  the  magnitude  and 
trends  of  the  two  sets  of  calculated  data  are  similar,  verifying  the  modest 
effect  of  the  changed  integration  limits  on  the  calculated  results.  The 
Figure  6  data  basically  depict  the  decrease  in  heterodyne  efficiency  with 
increasing  projected  target-radiance  solid-angular  subtense  (resulting  in  a 
smaller  radial  coherence  length) .  It  is  also  clear  from  the  longer-range  data 
of  Figure  6  (note  the  crossing  of  the  R  ■  5  km  and  R  *»  10  km  data)  that 
atmospheric  turbulence  eventually  becomes  the  dominant  factor  in  determining 
the  efficiency  after  the  solid  angle  is  reduced  sufficiently.  Thus  there  is 
an  optimum  range  for  maximum  detection  efficiency  depending  upon  the  target 
characteristics  and  turbulence  conditions. 

Tactical  Aircraft  Simulation 

To  provide  an  indication  of  the  heterodyne  detection  efficiencies  that 

might  be  expected  in  a  more  real-'  tic  laser-radar  detection  scenario, 

efficiency  calculations  were  perfo  ed  using  the  simulated,  tactical  aircraft, 

projected  cross-section  depicted  in  Figure  8.  The  Figure  8  cross-section 

distribution  assumes  a  spatially  constant  laser  bidirectional  reflectance 

-3  -1 

function  of  p'  -  5  x  10  sr  which  corresponds  to  the  sample  D  of  Reference 
11  (green  paint  on  metal)  at  a  backscatter  angle  of  70°.  This  reflectance 
function  yields  a  total  target  cross  section  of  48  cm  .  For  calculational 
purposes,  the  target  distribution  of  Figure  8  was  discretized  into  a  20  x  20 
grid  and  calculations  were  performed  as  a  function  of  range,  Figure  9,  for  the 
fully  illuminated  target  and  as  a  function  of  beam  size,  Figure  10,  at  4  km. 
The  transceiver  parameters  were  those  employed  for  many  of  the  previous 
calculations  and  were  chosen  to  be  representative  of  a  compact  installation. 

Two  simple  range  proportionalities  (R  and  R^)  are  depicted  in  Figure  9  to 
illustrate  the  fact  that  the  heterodyne-efficiency  range  dependence  decreases 
more  rapidly  than  R^  for  close  ranges,  increases  less  rapidly  than 
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Figure  8.  Simulated  tactical  aircraft  nose-on  projected  cross-section  The  laser  bidirectional-reflectance 
function  is  assumed  to  be  p'  =  5xl0'3  sr1.  The  superimposed  anulii  represent  the  assumed 
illuminating-beam  radii  used  for  heterodyne  efficiency  calculations. 


0  2  4  6  8  10 

Range  (km) 


Figure  9.  Calculated  heterodyne-detection  efficiency  as  a  function  of  range  for  the  simulated  tactical 

aircraft  target  depicted  in  Fig.  8.  Straight  lines  indicate  the  slopes  of  efficiencies  proportional 
to  R  and  R4  as  indicated.  The  fixed  parameters  of  the  calculation  are  A  =  P  =  2.5  cm, 
f g  =  fr  =  25  cm,  L0  =  5  cm,  X  =  10.6  fim,  and  C*  =  10‘14m'J/3. 
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Figure  10.  Calculated  heterodyne-detection  efficiency  as  a  function  of  illuminating-beam  radius  for  the 

simulated  tactical  aircraft  target  depicted  in  Fig.  8.  The  fixed  parameter  of  the  calculation  are 
A  =  P  =  2.5  cm,  f g  =  fr  =  25  cm,  L0  =  5  cm,  X  =  10.6  fim,  C2  =  10‘14  m’2/3,  and  R  =  4  km.  The 
illuminated-area  laser-radar  cross  section  is  plotted  and  the  product  of  the  calculated 
heterodyne  efficiency  and  laser-radar  cross  section  is  shown  as  a  relative  magnitude. 


R  for  long  ranges,  and  cannot  be  modeled  by  a  simple  power  relationship,  as 
one  might  expect  from  a  simple  coherence  area  argument.  The  rapid  decrease  in 
heterodyne  efficiency  for  close  ranges  implies  that  focusing  optics  may  be 
required  for  an  operational  large-aperture  system  to  maintain  an  acceptable 
signal-to-noise  ratio  because  the  signal  increase  with  decreasing  range  will 
be  proportional  to  only  (as  a  result  of  the  radar  range  equation).  The  net 
decrease  in  signal-to-noise  ratio  that  Figure  9  implies  for  close  ranges  may 
be  unacceptable  for  certain  applications.  Because  the  typical  dimensionality 
of  a  speckle  cell  decreases  inversely  with  the  range  (if  the  target  is  quasi- 
homogeneous)  one  might  expect  that  the  efficiency  should  decrease  inversely 
with  R^  because  the  efficiency  is  related  to  the  integral  of  the  coherence 
area  (Equation  (1));  however.  Figure  9  shows  that  this  is  not  the  case.  Also, 
Figure  9  shows  that  because  of  the  large  extent  of  the  target,  a  unity 
heterodyne  efficiency  is  not  achieved  at  any  range  because  turbulence  effects 
become  important  at  the  ranges  where  the  target  speckle  effects  become 
unimportant  (>  10  km).  Although  efficiencies  at  ranges  greater  than  10  km  are 
not  shown  in  Figure  9,  the  efficiency  decreases  at  ranges  much  greater  than 
10  km  because  of  atmospheric  turbulence. 
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Figure  10  demonstrates  that  the  heterodyne  eft*.ciency  can  be  increased  by 
decreasing  the  illuminated  target  area;  however,  the  converse  behavior  of  the 
laser-radar  cross  section,  also  plotted  in  Figure  10,  demonstrates  that  the 
net  gain  is  modest.  Because  the  actual  signal-to-noise  ratio  is  proportional 
to  the  product  of  the  heterodyne  efficiency  and  the  signal  power  (which  is 
proportional  to  the  laser-radar  cross  section),  the  product  of  the  heterodyne 
efficiency  and  laser-radar  cross  section  is  plotted  in  Figure  10  to  illustrate 
the  relative  actual  signal-to-noise  ratio.  The  Figure  10  results  indicate 
that  the  actual  signal-to-noise  ratio  does  not  decrease  beyond  a  beam 
divergence  ancle  of  0.25  mrad.  It  should  be  noted  that  the  constant  values  of 
heterodyne  efficiency  and  laser-radar  cross  section  indicated  for  the  smallest 
values  of  beam  radius,  result  from  the  finite  discretization  of  the  target 
cross  section  distribution  and  are  not  physically  meaningful. 

To  demonstrate  that  the  general  trends  illustrated  in  Figures  9  and  10 
are  not  unique  to  the  target  geometry  assumed  for  the  calculations, 
corresponding  results  are  shown  in  Figures  11  and  12  for  a  Schell-model  target 
that  would  correspond  to  a  uniform  diffuse  target  illuminated  with  a  Gaussian 


Figure  11.  Calculated  heterodyne-detection  efficiency  as  a  function  of  range  for  a  Schell-model  target 
with  a  1-m  amplitude  radius.  Straight  lines  indicate  the  slopes  of  efficiencies  proportional  to 
R  and  R4  as  indicated.  The  fixed  parameters  of  the  calculation  are  A  =  P  =  2.5  cm, 
f g  =  fr  =  25  cm,  L0  =  5  cm,  X  =  10.6  ^m,  and  Cj  =  10'14  m'2/3. 
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Schell-model  target.  The  fixed  parameters  of  the  calculation  are  A  =  P  =  2.5  cm,  f  =  ffr  ==  25 
cm,  L0  =  5  cm,  A  =  10.6  /xm,cj  =  10'*4  and  R  =  4  km.  The  laser-radar  cross  section  is 

proportional  to  the  Integrated  intensity  distribution,  which  is  proportional  to  the  square  root 
of  the  beam  radius,  which  is  also  plotted.  The  relative  magnitude  of  the  product  of  the 
heterodyne  efficiency  and  laser-radar  cross  section  is  also  illustrated. 

beam.  The  amplitude  radius  assumed  for  the  range-dependent  heterodyne- 
efficiency  calculation  shown  in  Figure  11  was  l  m.  Note  that  although  the 
magnitude  of  the  efficiency  values  shown  in  Figure  11  differ  from  Figure  9  (as 
one  would  expect),  the  functional  dependence  is  not  radically  different,  and 
similar  conclusions  may  be  inferred  regarding  the  necessity  of  focusing  optics 
and  the  ultimate  effects  of  atmospheric  turbulence  with  increasing  range. 
Similarly,  the  effects  of  beam  radius  shown  in  Figure  12  demonstrate  that  the 
actual  signal-to-noise  ratio  does  not  decrease  beyond  a  given  beam  radius  for 
the  Schell-model  target,  although  a  somewhat  larger  beam  divergence  is 
indicated  for  this  transition.  It  is  worth  noting  that  there  is  a  peak  in  the 
actual  signal-to-noise  ratio,  although  it  is  not  illustrated  in  Figures  10  and 
12,  because  the  signal,  perforce,  goes  to  zero  for  zero  beam  radius. 

Test  Panel  Calculations 

A.  test  panel  configuration  was  devised  as  a  simple  target  that  could  be 
used  to  validate  the  results  of  the  analyses  of  this  study  in  future 
experimental  tests.  Figure  13  depicts  the  suggested  panel  configuration 
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Figure  13.  Depiction  of  suggested  test-panel  configuration  consisting  of  either  four  10  cm  X  10  cm 
squares  located  at  the  corners  of  a  square  with  dimensions  up  to  2  m  X  2  m,  or  a  square 
of  dimensions  up  to  2  m  x  2  m.  The  panel  material  is  assumed  to  have  a  laser-radar  cross 
section  per  unit  area  of  10'2. 

consisting  of  either  small  square  panels,  located  at  the  corners  of  a  large 
square,  or  a  large  square  target  of  equivalent  dimensions.  Figure  14  shows 
the  results  of  heterodyne-efficiency  calculations  performed  as  a  function  of 
the  target  dimensions  assuming  the  transceiver  parameters  employed  for  the 
previous  calculations.  The  Figure  14  data  provide  an  Interesting  comparison 
between  the  extended  target  case  that  has  previously  been  analyzed,  and  the 
separated  target  case  (having  a  constant  laser-radar  cross  section  for  various 
separations)  that  obviously  gives  rise  to  a  loblng  structure  in  the  predicted 
heterodyne  efficiency  as  a  function  of  panel  separation. 

Several  conclusions  can  be  inferred  from  the  separated-targec  heterodyne- 
efficiency  curve  shown  In  Figure  14.  First,  because  the  calculated  negative 
heterodyne  efficiencies,  plotted  as  positive  values  in  Figure  14,  follow  a 
continuous  curve  with  the  calculated  positive  values,  it  is  reasonable  to 
assume  that  they  are,  indeed,  meaningful  efficiency  values.  Thus,  the 
physical  interpretation  of  negative  efficiency  values  (discussed  in  Section  5) 
as  representing  a  quadrature  component  is  validated  to  the  degree  that  no 
other  physical  interpretation  seems  reasonable. 
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Figure  14.  Calculated  heterodyne-detection  efficiency  as  a  function  of  target  extent  for  the  two  test-panel 
configurations  shown  in  Fig.  13.  The  fixed  parameters  of  the  calculation  are  A  *  P  ■  2.5  cm. 
fj>  =fr  =  25  cm,  L0  =  5  cin,  A  =  10.6  ^m,  =  10‘14  m'2/3,  and  R  -  1  km.  The  data  plotted  as 

( * )  are  calculated  values  of  the  heterodyne  efficiency  which  are  negative,  but  have  been  plot¬ 
ted  as  positive  quantities. 


Second,  the  lobing  pattern,  observed  in  Figure  14,  is  redolent  of  a 
fringe  visibility  pattern,  such  as  one  might  obtain  in  a  coherence 
experiment.  Indeed,  coherence  effects  are  responsible  for  the  degradation  in 
heterodyne  efficiency  and,  because  of  the  correspondence  of  speckle  effects 
with  limited  spatial  coherence®,  it  is  not  surprising  that  separated  objects 
should  give  rise  to  a  fringe-like  visibility  pattern  at  the  focal  plane  that 
varies  with  target  separation  in  an  analogous  fashion  to  the  fringe  visibility 
technique  used  to  measure  stellar  diameters^.  This  correspondence  can  be 
demonstrated  by  considering  the  inverse  problem  of  measuring  the  receiver 
aperture  (now  interpreted  as  an  incoherent  disk  source)  by  measuring  the 
fringe  visibility  at  the  target  plane.  The  first  null  in  Figure  14  is 
observed  for  a  panel  separation  of  1.2  ra,  which  would  imply  a  disk  size  of 
approximately  1.1  cm  for  the  Inverse  problem  If  the  fringes  were  obtained  from 
two  point  sources.  The  fact  that  the  panels  are  finite  and  arranged  in  a 
square  fashion,  coupled  with  other  optical  factors,  accounts  for  disagreement 
between  the  prediction  and  the  actual  aperture  diameter  of  5  cm.  Calculations 
of  the  heterodyne  efficiency  for  two  small  panels  (analogous  to  the  pinholes 
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tn  a  Young  Interference*  experiment)  yields  a  null  for  a  separation  of  0.25  ra, 
which  corresponds  to  an  aperture  diameter  of  5.2  cm  for  the  inverse  problem, 
which  Is  in  good  agreement  with  the  actvial  aperture  diameter. 

Finally,  it  Is  noteworthy  that  the  loblng  pattern  shown  in  Figure  14 
would  be  replicated  In  the  measured  s Ignal-to-nolse  ratio  of  a  laser-radar 
system,  as  well  as  the  heterodyne  efficiency,  because  the  laser-radar  cross 
section  of  the  target  does  noi  vary  with  separation.  Thus,  the  suggested  test 
panel  configuration  provides  a  relatively  simple  means  of  testing  some  of  the 
key  predictions  of  the  analysis,  namely  (a)  coherence  effects  and  (b) 
quadrature  components. 


32 


7.  GENERALI7.R0  laser-radar  heterodyne-efficiency 
COMPUTER  CODE 


The  computer  code  LADAK  that  calculates  the  laser-radar  heterodyne 
efficiency  for  either  a  Schell-model  target  or  general  target  cross-section 
distributions  (assuming  either  a  circular  aperture  or  Cassegrain  receiver 
geometry)  Is  listed  below.  The  code  listing  is  preceded  by  a  brief  exposition 
of  the  input  parameters  required  for  its  operation.  All  Input  parameters  are 
in  31  units.  The  code  assumes  a  laser  wavelength  ot  10.6  pm,  which 
corresponds  to  the  peak  in  the  gain  curve  of  the  P  branch  [P  (?.0)]  or  the 
(00°1)  -  (10°0)  band  of  the  CO2  laser.  Alternative  laser  transitions  can  be 
studied  using  the  LADAR  code  by  editing  the  value  assumed  for  the  variable 
LAMBDA  in  an  appropriate  manner. 


LADAR  Input  Parameters 


SIuPHI 


-  Local-oscillator  aperture  radius. 

-  Entrance-pupil  inner  radius.  A  non-zero  value  indicates 
Cassegrain  geometry. 

-  Entrance-pupil  outer  radius. 

-  Receiving  lens  focal  length. 

-  Local-oscillator  lens  focal  length. 

-  Local-oscillator  Gaussian  amplitude  radius. 

-  Range  to  target, 

-  C2 

n 

-  Phase  variance  <j.  . 

<P 

-  Glint-poiut  Gaussian  radius  of  curvature. 

-  Source  amplitude  radius  L^.  A  non-zero  value  indicates  a 
Schell-model  approximation. 

-  Average  phase-front  radius  of  curvature  L<$>  for  the  Schell- 
model  approximation. 

-  pu  for  the  Schell-model  approximation. 


f 
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If  LA  **  0,  the  following  are  required. 

M1,N2  -  Number  of  equally  spaced  integration  points  used  for  the 

numerical  integration  over  the  two  components  of  Z,  . 

The  following  parameters  describe  the  target  laser  scattering  distribution. 


CPIS,  ZP2S 

ZP1E.ZP2E 

NP1.NP2 

D2S 


-  Lower  bounds  of  the  scattering  distribution  for  the  two 
Cartesian  components  of  £  . 

-  Upper  bounds  of  the  scattering  distribution  for  the  two 
Cartesian  components  of  (  . 

-  Number  of  equally  spaced  integration  points  used  for  the 
numerical  integration  over  the  two  components  of  £  . 

-  NP1  x  NP2  matrix  of  scattering  distribution  values. 

at 


matrix  of  scattering  distribution 

D2S  (IP1.1P2)  is  the  value  of  — — 

dC 


ZP1S  + 


ZP2S  + 


(2*IP1  -  1)  (ZPIE  -  ZPIS) 
2  NP1  1 

(2*IP2  ~  1)  (ZP2E  -  ZP2S) 
2  NP2 
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«F  TX?n 

REAL  LAMBDA, K, LO, LA, LPHI 
COMMON  /WAVE/  K 

COMMON  /REC/  A/ P  IN/ PCUT, FP,FL,LO 

common  /tar/  r,cn?,la, lph,rhou 

COMMON  /CROSS/  N«>1,  Z  PIS,  ZP  IE  ,NP?,  7.  »2 S,  ZP2f  ,02$  !  21,  21) 

OATA  PI  /3. 14159265/ 

LAMBDA-10. 6E-6 
K-2. *P!/LANBDA 
WRITE! 6# 1 ) 

1  FORMAT!*  ENTER  RECEIVER  PARAMETERS;  A. RT N# POUT, F R* F L, LO • ) 
»EADI5,)  A,PIN,0OUT,FP.»FL.»LO 

WRITE!  6,  2 ) 

2  FORMAT!*  ENTER  TARGET  PARAMETERS?  p»  CM  2,  S  IGPHI,  RADG*  ) 
READC5,)  R,CN2#SIGPHI, RACG 

WRT  TE! 6, 6 ) 

6  FORMAT!*  ENTER  TARGET  PARAMETERS?  L A, l RH I, RHOU* ) 

LA-0.0 

RPAD!5,»  LA, LPHI. RHOU 
C  CALCULATE  THE  COHERENT  EFFICIENCY. 

CALL  ECOH(EFEC) 

IF!L A.EQ.O.O)  WRITE!6,11?  EFFC 
11  FORMAT!*  COHERENT  EFFICIENCY  «'»Gl?.4? 

Ip!  L  A.NE  •  0  •  0  '  WR  ITE  ( 6. 1 5  I  CFFC 
IE(L A.NE.O.O)  STOP 
A»-A*FR/FL 
RMIN-AMIN! AP.PIN) 

RM0UT-AMIN!AP,P01T) 

WRI T E ! 6,  12  )  RMIN.PMOLIT 

1?  FORMAT!'  RMIN  -  ' G 1 2  .4, /,  •  PMPIJT  »*,G12.4> 

C  READ  INTEGRATION  PARAMETERS. 

WPIT  F( 6*  3) 

3  FORMAT!'  ENTER  1  INTEGRATION  PARAMETERS;  N1,N2«) 

READ!  5,)  N1,N2 

C  READ  SCATTERING  CI5TP1BLTI0N . 

WR I TF ( 6, 4 ) 

4  FORMAT!"  ENTER  Z"t  INTEGRATION  P AR  AMET  F&S*  Z*1S » 7P1E, NP1  •  , 
1  /, 3 4X, " 7P2S, ZP 2F, NP2 «  ) 

READ  !2, )  7P15,ZP1E*NP1 

n7Pl»!7PlE-ZPlS)/NPl 

READ  !2* )  ZP2S.ZP2E »NP2 

0ZP2«!ZP2F-7P2S»/NP2 

WRITE! 6, 5)  . 

5  FORMAT!*  ENTER  SCATTERING  DISTRIBUTION?  D2S*  ) 

DO  60  IP2-ls»NP2 

RFADI2, )  !D2S!IP1,IP?),IP1-1,NP1) 

60  CONTINUE 

C  INTEGRATE  THE  SCATTERING  DISTRIBUTION.! 

SIGTAR-0,0 

Zn-ZPlS+DZPl/2. 

00  200  IP1-1.NP1 
7P2-7P2S+DZP2/2. 
nn  150  IP2-1.NP2 
STGTAR-SIGTAR+D2SIIP1, IP2) 

150  ZP2-ZP2+DZP2 

?°°  5liT56^?IfiRTDZPlTDZP2 

WRITE!6, 13)  STGTAR 

13  FORMAT!  *  TOTAL  TARGET  INCOHERENT  I.RCS  •*»G12«4) 

C  CALCULATE  THE  INCOHERENT  EFFICIENCY. 
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CALL  E INC0M<N1#N2*PFFI) 
EEFI-EFFI/SrGTAR 


WRITF(6,141  fffi 

FOR  MAT  { *  INC  OH  ERF  NT  EFFICIENCY  «'»G12.'4l 
TfiRMi-<XEXP<S  IGPFI»S  IC-PHII-l.)*SIOTAP 
TERM?-PI*RAD-I  *RAOG 
E°S»0.0 

IP!  TFRM2.NE  .0.0)  E  PS- E  P  S  ♦  Ef  FC  /  ( 1  .  *T  F  RM  1/’ T  E  RM2 } 
Ic<  TFRMl.NE.C.O)  fPS-EPS«EFFim.  +  Tf  rM^/tfRMI) 

IF ( f PS  »NE • 0  «  0 )  WPITF(6»15)  EPS 

FORMAT  < 1  MAXIMUM  HETERODYNE  EFFICIENCY  -',G12.4) 

STOP 

£N0 


SUBROUTINE  ECOM{ Ef FC ) 

REAL  QC(3) 

REAL  LP,K.LO,LA,LPHl 
real  lap,lphip 
COMMON  /WAVE/  K 

COMMON  /REC /  A, P  IN, POUT, f R, FL,LO 
COMMGN  /TAR/  R , C N2, L A# L P H I , PHOO 
COMPLEX  AL°HA,BETA 

Complex  xcexp 

COMMON  /FINT/  Al PH A, RMOP , R M 
COMPLEX  IC 
ev^ernal  F 

OATa  PI,IC  /3.1415<J2fc5,(C.,l.>/ 

RHOO-1  . E  6 

IF( CN2.NE.0.0)  RHQO- <3.*1. 455*K*K *» *CN 2 /8 . >*♦(-. 6) 
RHOP-RHOO 

C  CHFCK  FOP  SCHELL-MODEL  CALCULATION. 


I*(L A.EO.O.O)  GO  TC  10 
ALPHA-CMPLX(2./K/LA» 1./LPHII 


ALPHAI-AIMAG(ALPHA) 

Z2-2.*(2.*R/K/LA  1**2 

RHUP  -S  QR  T  ( (Z2*(RHG0**(-2  J  +  RHQU** ( -2 ) )♦ C ABS < 1 . ♦ XC*AIPM A+R )**2 ) / 

1  (  RHOO**  <-?)*<  1.-2.  *R*Al  PHAI  +  CARSI  ALpO*R  )**2)  + 

2  RHO 0** ( -2  )*  IP*ALPHAI*,75*22*RHO0*M-2  M  ♦ 

3  RHOU**(  — 2)*(1.^PHC0’*Y(-2)*Z2)H 

L4P-LA*$QRT<  Z2  + ( RHCO ♦* (- 2 ) ♦RHOU» *  I -2 ))  ♦CABS! 1. ♦ IC* AL PHA*R )**2) 
Lt>HIP»R*(Z2*(RH00**  (-2  )  ♦  RHOU  **  1-2  >  )  ♦  CA  B$  ( 1  .♦  IC*  AlPH  A*R  )  **2 )  / 

1  (ALPHAI*R-CABS(ALPHA*P)**2-Z2*(  , 5* R HO 0* * ( -2 ) ♦RHOU** (-2 )»> 

10  CONTINUE 

AP«A*FR/FL 

Z«K/2.*(1./FL-1. /FP1 

ZP-Z*(Fk/K)**2 

LP-LO*FP/FL 

BETA -CMR LX (l P*LP,-4.*ZP  J  /  ( L P**4  + 1  & . *7P ** ? ) 

TF(L A.EQ.O.O)  AL PH A- BET A +C M PLX < P HOP **< -2 ) ,  K/2./R ) 

|E( LA.NE .0.0)  ALPHA«PETA*CNPLX(PH3p**(-2M-X./LAP,K/2. /IPHIP) 
ScTAP-REALlBETA) 

PMIN-AMIN(AP,PIN) 

RMOUT«AMIN(AP,POUT ) 

DO  100  1-1,3 

0C(  n-o. 

PM1-RMIN 
RM2 - RMOUT 

C  CHECK  FOR  CASSEGRAIN  RFCFIVEP  GEOMETRY. 

IF( RMIN.EO.0.0  .ANC.  I.GT.ll  GO  TO  100 
IF(I.EO.l)  RMl-RMOtT 
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IF(  I  . E 0 , 3 )  RK,2-Rhlfs 
C  CALCULATE  THE  COHERENT  INTEGRAL. 

XO  E 1  -1  .-XC  EXP { -ALPHAS*  1«RM1)-*CEXP(-C  QNJG(A(.PHA)*RM2*RM2) 
XCE2-XCEXP  (- ALPHAS PM1*PM1-CON.JO(  iLPHAT*RH2*RM2  + 

1  2.*PH1*RH2/RH0P/RHCP) 

or  (I  )-(XC61*XCE2*tIC<2.*PMl*RM2/RH0P/RH0P|  )  /  A* /CABS  (  ALPHA)**2 
PM-R*2 

CML  SIMP<F*0.#RMl#0.01*10#SP#SF#Nt, IER) 

I«=(PP1.EQ.RN2)  S  F-2 • C*S  F 

q2H1*QC(I}-RM2*SF/2./  (CABS(ALPHA»»RH0P»A#2 

T*<PMl.EQ.RM2)  GC  TO  50 

PM-PM1 

CALL  SlMp<  e,0.#RM2#0.01#10#SP*SF#NT,  IER1 

.0  SSHi:8SIK;S9tKaii»tfUna!Mn!SM»iHari**fi 

100  CONTINUE 

£cFC«?.*8eTAP/t®0UT*«2-PIN*A2)*(QCm-2.*QC<2»*QCC3i  1/ 

1  (XEXP(-2.*B£TAR*RMIN**2)“X€XP(-2.I*BETAR*R MOUT+*2| ) 

RFTtIRN 

END 

C  SUBROUTINE  E INCQH ( Nl#N2#  EE  FI )  ' 

PCAL  QI(3) 

COMMON  /RFC/  A#  P  IN#P0UT.#f  R#F  L#L0 
rOMMCN  / TAR/  R,CH2« LA*  L«»H#RHnu 

COMMON  /CROSS/  NP1* ZP1S* ZP1E*NP2# ZP2S. ZP?F >D2SI 21, 21 ) 

COMPLEX  AL  PH  A#  BETA 
PEAL  JINC 
COMPLEX  1C 

OATA  PI,IC  /3.L415<!265,  (C.#l.)/ 


COMPLEX  !( 
OATA  PI,I( 
RHnO«l.E6 


IP(CN2.NE.0.0)  RHQQ- < 3  *  *  1  • 455*K*K*R  +  CM  ? / R . )**(-• 6) 
PHQP  »R  HOO 
AP- A*FR/FL 

7-K/2.M1./FL-1./FR) 

ZP-  Z  * ( ER/K )  +*2 
L°-LO*FR/FL 

BcTA»CMPLX(LP*LP>-A.*ZP)/(  LP**A«-  l*»,l*ZP**2  J 
ALPHA«BETA«-CMPLX<RH0P**<-2T#K/?./R) 
BETAR-REAL(BETA) 

BPTAT-AIMAG1BETA) 

BETAT-BETAR  +  2./RH0P/RH0P 
RMIN-AMIN( AP*PIN) 

RMOUT-AMINt  AP,POUT) 

C0NST-2.*P/KA>BETAT  +  1.0 
00  5  CO  1-1# 3 
QIC  I )" 0 • 

Pn- R  NOUT 

CHECK  FOP  CASSEGRAIN  RECEIVER  GEOMETRY « 

IF(R*‘1N.EQ.0.0  .ANO.  I.GT.l)  GO  TO  500 
IF(I.E0.2)  PP-{ RMIN+RM0UTI/2. 

IF(  I  .E0.3  )  RP-RMIN 


■BFTAR+2./RH0P/RH0P 
iMIN(  AP#PIN) 


GEOMETRY. 

,1)  GO  TO  500 


IECI.E0.2)  PP-(RMIN+RM{HJT)/2. 

IF(  I  .E0.3  )  RP-RMIN 
PM-PP 

Ic(  I.EQ.2)  RM-SQRT  ((PNIN*A>  2+ R  MOOT  **2 )  /  2.  I 
CALCULATE  THE  INCOHERENT  INTEGRAL. 

On-  2,  ♦RP/Nl 

Z1S--PP 

0Z2-2.*RP/N2 
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-  i«ui 


c  — 


7?S--RP 

P?pi«(zpie-ZP1S)/NP1 
GZP?-( ZP2F-ZP2S )/NP2 
FPS-0.0 
Z1-Z1S+DZ1/2. 

DO  4CQ  11*1, N1 
72-Z2S+OZ2/2. 

DO  300  12-1, N2 

IE<Zl*Zi*Z2*Z2.GT.PP*RP)  GO  TG  300 
Z°l-ZPlS+D7Pl/20 
on  200  I  P1-1.#NP1 
Zp2-ZP2S*DZP2/2» 

OG  ICO  IP2-1,NP2 

IF(D2S(IP1*IP2) .EG. 0.01  €0  TO  100 

U®-S  GRT I ( Zl*C0N5T+  7  Pll**2* (Z2*C0N? T ♦Tp21**2)*K/R 

A«2.*SCRT(RM*PM-Zl*21-Z2*Z2i 

IF(UP.LT.BETAT*A>  GC  TO  «o 

FT- 2  .^P I* A* A* J INC (UP* A i 

GG  TO  90 

FT-?.*pI*<  XEXP(-UP*UP/?«/R£TAT)/RFTAT-XFXP(-A*A/2.*( 

1  BETAT*UP*UP/2. )  )/ < BfTAT*UP*Up/? . ) ) 

£°S»Jp^£xEXP<-2.*BETAR*(Z1*Z1*7?*Z?)  M'FT+D2S(Ipl»IP2> 

Z®2-  ZP24-0ZP2 
7H-ZPU0ZP1 
Z7-72+DZ2 
Z1-Z1+DZ1 

QI(I)-El>S*DZl*DZ2*DZPl*0ZP2 

CONTINUE 

EEFI-BETAR/PI/Pl/(POUT+*?-PIN**?)*<QTt  1» -2  •* 01  (2  ) ♦01  «  3)  »  / 
1  (XEXP  <-2.*BETAR*RMIN**2)-XEXpt-2.*BFTA»>*RMnUT*A2n 
RETURN 
END 


FUNCTION  F(Y) 

COMPLEX  ALpH/,XCEXP 
COMMON  /TINT/  ALPHA, PC, PM 

F«XCEXP(-ALPHA*YAY-CCNJG (AlPHA1*RM*RM+2.*Y*RMyRC/RC) 

IF(F.EO.O.O)  RETURN 

F«F*£I1(2«*Y*RM/PC/RC) 

RETURN 

END 


COMPLEX  FUNCTION  XCEXPOO 
COM® IE X  M 
XCEXP-  ( 0  •# 0 •  I 
IMREAL(W) .LT.-650.)  RETURN 
XCEXP-CEXP(W) 
rcturn 
END 


FUNCTION  XEXP(X) 
•'6XP-0.0 
IF<  X.LT.-650.) 
XtXP-EXP  (XJ 
RPTurn 


RETURN 


PEAL  FUNCTION  JINC<X) 
IC(X.GT.3.0)  GO  TO  ICO 
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T*( X  /3.0)*  +  2 

JINC-0.5+T*  (-0.5 6249985  + T+( 0.21093573+  T*  (-0. 0305*? 89+ 

1  T*  (0.00443319+T*(-0.00031761+T*Q.00001109)>  ))) 

rei  urn 

100  T-3.0/X 

F1*0«797684'56+T*(0.CC0CC156+T*(0,01659  66  7+T*(0«  000 17105+ 

1  T*(-0 .002 44511  + T*(C. 0011 3653-T*0.  0002  Of)  33  }  II )  ) 

THE TAl«v-2. 3  5619449  +  T*(C.l 249961  ?♦ T* (0 .00005650+T* < 

1  -0, 006  :>7  97  4+T*(0.CCO74  3  4«+T*(0, 0007  98  ?4 -T*0. 00029166  >  )  )  )) 
JTNC-F  1*C0S  (THETAl  WX  +  M  1.  5) 

RPTUPN 

ENn 

FUNCTinN  fjoTxT  **  *  —  — 

C  COMPUTE  EXP(— X )  ♦  10. 

IF( X.GT.3.75)  GO  TC  100 
T-(X/3.75>**2 

ETO-  l.+T*(  3. 5156 224  + T*(  3. 0  699424  + T*(  1.  206749 ?+T*( 

1  ,2  659732  +  T* (.0360768+T*  ,0045813)1  ) ) ) 

EI0-FXP(-X)*fI0 

RFTU&N 

100  T-3.75/X 

ET0».39894228+T*{«01328592+T*(.002?5319+T* C-. 00157565+ 

1  T*(.0 09 16281 +T *(-.02057706 +T*(.Q263 55 37+T*( -.01647633+ 

2  T*. 00392377))))))) 

P!0-EI0/S0RT(X) 

RETURN 

END 

c  funct7on”eii77) 

C  COMPUTE  EXP(-X)  +  II. 

TFJX.GT.3.75)  GO  TC  100 
T-(X/3.75)**2 

Ell-C. 5+T*( .0769C594+T*(  .51498869  +  T*(?1508  4934  +  T*( 

1  , 02 65 87 33 +T*(. 00301 53 2+T*. 00032411) )) )> 
fc  I  1  *  *"  I  1*  X 
EI1-EXP(-X)*EI1 
RETURN 

100  T-3.75/X 

ET1-  . 39894 2 2 8+T* (-. 03988C24  + T*( -. 00362 01 8  +  T*(. 00163 801  + 

1  T* (-.01031555 +T *(.02282 96 7+T*(-.0?895312+T*t 

2  .01787654+TJ  (-.  C04  2C0  59  )))))))  ) 

EI1*EI1/S0RT (X) 

RETURN 

E*'0 

"SUBROUTINE  $If«P(F*  A*fi*0EL*  IMAX*Sn,S  JnJTerT" 


EIMPSON 

N-0 


INTEGRATION 


INE  • 


BA«B-A 

IF  (  8A ) 20* 19  »  20 

IpR  « 1 

RETURN 

I  e  ( DEL ) 22*  22*23 

IER-  2 

RETURN 

IF(lNAX- 1)24*2 4> 25 
IEP  ■  3 
RETURN 
X-BA/2.  +  A 
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